"Quasi 2-D" Spin Distributions in II- VI Magnetic Semiconductor Heterostructures: 

Clustering and Dimensionality 
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Spin clustering in diluted magnetic semiconductors (DMS) arises from antiferromagnetic ex- 
change between neighboring magnetic cations and is a strong function of reduced dimensionality. 
Epitaxially-grown single monolayers and abrupt interfaces of DMS are, however, never perfectly 
two-dimensional (2D) due to the unavoidable inter-monolayer mixing of atoms during growth. Thus 
the magnetization of DMS heterostructures, which is strongly modified by spin clustering, is inter- 
mediate between that of 2D and 3D spin distributions. We present an exact calculation of spin 
clustering applicable to arbitrary distributions of magnetic spins in the growth direction. The re- 
sults reveal a surprising insensitivity of the magnetization to the form of the intermixing profile, 
and identify important limits on the maximum possible magnetization. High-field optical studies of 
heterostructures containing "quasi-2D" spin distributions are compared with calculation. 

PACS numbers; 75.50.Pp, 78.55.Et, 71.35.Ji, 75.50.Cn 
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Spin clustering is ubiquitous in II- VI diluted mag- 
netic semiconductors (DMS), resulting in reduced effec- 
tive magnetizations at low magnetic fields and magne- 
tization steps at high fields.El Clear predictions can be 
made for the number and type of spin clusters in 3D sys- 
tems {e.g., bulk Cdi-ajMn^jSe), where the distribution of 
magnetic Mn^+ cations is random and isotropic. With 
the advent of molecular beam epitaxy (MBE) and other 
techniques for monolayer-by-monolayer growth of DMS 
heterostructures, spin clustering in systems of reduced 
dimensionality hajS enjoyed much recent interest. It is 
well establishedo^El that spin clustering (arising mainly 
from an antiferromagnetic exchange between neighbor- 
ing magnetic cations) should be greatly reduced in two- 
dimensional systems such as abrupt interfaces or discrete 
monolayer planes, leading to enhanced paramagnetism. 
However, experiments show that perfect 2D interfaces 
and monolayers are never realized due to the inevitable 
inter-monolayer mixing of atoms during MBE growth, 
which smears the magnetic cations over several mono- 
layers. Common mechanisms include segregation (mix- 
ing between the monolayer being grown and the underly- 
ing monolayer) which leads to roughly exponential mag- 
netic profiles, and diffusion (which can arise from, e.g., 
high growth temperatures or annealing) which leads to 
gaussian profiles. Hence, real DMS heterostructures are 
more accurately said to contain "quasi-2D" distributions 
of spins, with a corresponding magnetization and degree 
of spin clustering somewhere between that of bulk (3D) 
and planar (2D) spin distributions. 

The local, planar magnetic concentration in these 
quasi-2D spin distributions varies significantly from 



monolayer to monolayer, strongly affecting the probabil- 
ity of forming spin clusters (which themselves may span 
many monolayers). It is desirable to quantitatively pre- 
dict the degree of this spin clustering in a given DMS het- 
erostructure so that accurate comparisons can be made 
with real data. In this paper we present exact expressions 
for determining the number and type of spin clusters (sin- 
gles, pairs, open- and closed triples) for arbitrary distri- 
butions of magnetic spins in the common (100) growth 
direction. The results reveal a rather surprising insensi- 
tivity of the computed magnetization to the form of the 
intermixing profile (exponential/gaussian), and highlight 
important limits on the maximum possible magnetiza- 
tion using MBE techniques. High-field photolumines- 
cence (PL) and reflectivity studies of DMS superlattices 
and quantum wells containing quasi-2D magnetic planes 
are compared with the analytic results. 

Spin clustering in DMS derives predominantly from the 
strong antiferromagnetic d-d exchange between nearest- 
neighbor (NN) magnetic cati«ns ( Jata? ^ — lOif). As 
outlined in the work of Shapirael and otherg3, single Mn^"*" 
cations with no magnetic NNs are S = ^ paramagnets, 
with Brillouin-like magnetization. Two NN Mn^+ cations 
form an antiferromagnetically-locked pair with zero spin 
at low magnetic fields, and step-like magnetization at 
high fields and low temperatures. Three Mn^+ spins can 
form a closed or open triple with net spin St = ^ and 
St — \ (respectively) at low fields, and a unique set of 
magnetization steps at high fields Jj Spins in higher order 
clusters are usually treated empirically and oftjen exhibit 
a linear susceptibility at high magnetic fields E 
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FIG. 1. Measured energy shift (abs. value) of the exciton 
PL to 60 T at 350mK in three quantum weUs containing the 
same number of Mn^"'" spins, but in different quasi-2D distri- 
butions, showing increased spin clustering with increasing 2D 
concentration. 

The magnetization of monolayer planes of Mn^+ spins 
is a significant challenge for conventional magnetometry. 
Alternatively, the magnetization from DMS heterostruc- 
tures may be inferred from their giant magneto-optical 
properties. The Jsp-d exchange interaction between elec- 
trons/holes and local Mn^+ moments generates giant 
exciton spin-splittings that are proportional to the lo- 
cal magnetization within the exciton wavefunction. Us- 
ing the giant spin-splitting of confined excitons to probe 
the jjnagnetization. within a quantum well, the studies of 
GajH, Grieshaberl, and of OssauQ clearly established i) 
an enhanced paramagnetic response in very thin layers of 
magnetic semiconductor, and ii) that "ideal" magnetic- 
nonmagnetic semiconductor interfaces are smeared out 
due to segregation of Mn^+ during growth. A clear ex- 
ample of both effects can be seen in the high-field PL 
data of Figure 1. Here, we measure the giant energy 
shift (oc magnetization) of the band-edge exciton PL to 
60 Tesla in three quantum wells, each containing the 
same total number of Mn^+ spins, but in very different 
distributions. The samples are 12QA ZnSe/Zn.8Cd.2Se 
single quantum wells into which the magnetic semicon- 
ductor MnSe has been incorporated in the form of 'dig- 
ital' planes of 25%, 50%, and 100% monolayer cov- 
erage (12, 6, and 3 planes, respectively). The sam- 
ples and tie experimental method have been described 
elsewhere.Q Evidence of decreased spin clustering with 
decreasing planar concentration is clear in the low-field 
magnetization {H < 8T) , which is largely due to isolated 
Mn^+ spins and is fit to a modified Brillouin function, 
Esat = B;^/2[5^BH/kB{T + Tq)]. Esat is the saturation 
splitting and Tq is an empirical parajneter that accounts 
for long-range Mn-Mn interactions.^ As shown, as the 
planar Mn^+ density increases, Esat decreases while Tq 
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FIG. 2. Probability of clustering for a spin in the n*'' (100) 
monolayer of an arbitrary Mn^"'" distribution. Multiplying by 
Xn gives the number of clustered spins. Cations in the n + 1*'' , 
n"*, and n — l"* monolayer are labeled by white, shaded, and 
black dots respectively. 

grows, consistent with the expectation that increasing 
the planar density results in fewer single (un-clustered) 
Mn^+ spins and more long-range correlations between 
Mn^+ spins. Further, with increasing planar density, the 
high-field susceptibility evolves from magnetization steps 
(from Mn-Mn pairs), to the linear susceptibility common 
in large, highly correlated spin clusters. The presence of 
low-field paramagnetism in the 100% MnSe planes sug- 
gest, however, the quasi-2D nature of these spin distri- 
butions: a full (perfect) 2D magnetic monolayer would 
behave as one infinite cluster and contain no paramag- 
netic spins whatsoever. 

The probability of forming a particular spin cluster is 
essentially determined by the number of (possibly mag- 
netic) NNs bordering the cluster, and the number of ways 
the cluster can form. In the zincblende crystal structure 
we consider, cations form an fee lattice, with twelve NNs 
per cation. Thus the probability of a Mn^+ being single 
or paired in bulk crystals (3D) is P^^ = (1 — and 
Pp^ = 12a;(l — x)^^ respectively, where x is the Mn^+ 
concentration. For perfect 2D monolayers grown in the 
(100) direction, the cations form a 2D square lattice 
with only four possible magnetic NNs per cation, so that 
p2£i = (1 - x)^ and P2£i ^ 42.(1 _ 2:)6. in a real sys- 
tem, the effects of diffusion and/or segregation intermix 
adjacent monolayers, so that a "perfect" 2D plane of 
DMS is smeared over several monolayers, with the n*'' 
(100) monolayer having a 2D magnetic concentration Xn 
(assumed to be random within the plane). Clustering 
within these. auasi-2D spin distributions can be mojieled 
numericallyH'Q or through analytic approximationsH, but 
an exact expression has been lacking. Figure 2 shows 
the exact probabilities of a Mn^+ spin in the n*'' mono- 
layer belonging to a single, pair, closed- and open triple. 
Diagrams show the clusters under consideration - e.g., 
there are three different types of Mn^+ pairs (with the 
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FIG. 3. Schematic of the Mn^^ concentration, including 
segregation effects, in a 10 monolayer wide quantum well 
with a::A/n=30% magnetic barriers, b-e) The calculated num- 
ber of Mn^^ cations (per thousand sites) in singles, pairs, 
open/closed triples, and higher order spin clusters. 



paired spin in the n — 1*'', n*'', or n + 1*^ monolayer), 
each four-fold degenerate. There are four types of closed 
triples (for a total of 24), and 126 total configurations 
for open triples. (We do not attempt the 1900 configura- 
tions of spin quartets that have been recently identified in 
the bulkllj, nor do we consider the much weaker distant- 
neighbor couplings between Mn^+ moments.EJ) 

This algorithm allows for an exact calculation of spin 
clusters in a heterostructure with an arbitrary distribu- 
tion of Mn^+ in the (100) direction. An example of its 
utility is shown in Fig. 3, where we compute the num- 
ber of Mn^+ in singles, pairs, triples, and higher-order 
clusters in a 10 monolayer (10 ml) wide nonmagnetic 
quantum well with magnetic {xmu = 30%) barriers. We 
assume full segregation of atoms during growth, giving 
e""/^ and 1 - e""/^ Mn'^+ profiles (A - l.Uml) at the 
first and second interfaces (Fig 3a). Figs. 3b-e show the 
type and number of clusters in each monolayer. Although 
the Mn^+ density is comparatively small near the center 
of the quantum well, the paramagnetic contribution from 
single spins to the optically-measured magnetization can 
be significant, as it depends on their overlap with the ex- 
citon wavefunction as shown. The density of triples and 
of pairs is clearly peaked in the quasi-2D interfaces. In 
the barriers, the spin distribution is bulk-like and the vast 
majority of spins are bound up in higher order clusters 
that contribute little to the low-field magnetization. 

Modeling quasi-2D spin distributions leads to some 
rather unexpected results. In particular, it is clear that 
magnetization studies alone will be of limited use in dis- 
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FIG. 4. a) Calculated average magnetization (per spin) 
for the Mn'^"'' spin profiles shown. Despite differing inter- 
mixing profiles, the magnetization is often indistinguishable, 
b) Calculated number of isolated Mn^"*' spins (per thousand 
sites) assuming a 8% bulk spin distribution is redistributed 
as digital planes (see text). Only if planes are spaced every 
other monolayer is there an enhancement, although segrega- 
tion greatly reduces the effect. 



tinguishing the exact form of the intermixing profile. Fig. 
4a shows the calculated magnetization for four differ- 
ent profiles of an initially 2D monolayer containing 20% 
Mn^+, where we include the magnetization from singles, 
pairs, triples, and higher order clusters following Refs. 
1 and 6. Though unrealistic, the first two profiles - 
a perfect 2D plane with xmu = 20% and two adjacent 
planes with xmu = 10% - illustrate an important point: 
clustering often "conspires" to equalize low-field magne- 
tizations. Although the single monolayer contains 5% 
fewer single Mn^"*" spins, it contains over a third more 
open triples and higher-order clusters, which act to equal- 
ize the deficit. Only at the first magnetization step are 
the profiles distinguishable, as the single monolayer con- 
tains fewer Mn-Mn pairs. The last two profiles represent 
the exponential and gaussian profiles roughly expected 
from segregation and diffusion, respectively, with decay 
length and half-width equal to 1 ml. Again, the calcu- 
lated magnetizations are nearly identical (although larger 
than for the first two profiles). Thus, magnetization mea- 
surements alone cannot distinguish the form of the spin 
profile. However, assuming a particular form, the mag- 
netization does depend sensitively on the segregation (or 
diffusion) length, which can then be used to fit an inter- 
mixing lengthscale as demonstrated below. 

The model we present can also identify configurations 
for realizing the maximum possible magnetization per 
unit volume in MBE-grown structures. One motivation 
for growing 'digital' alloys is to exploit the reduced clus- 
tering of 2D planes to achieve enhanced magnetizations 
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beyond those possible with bulk, 3D distributions.l3't3 In 
bulk DMS, the maximum paramagnetic response is ob- 
tained with XMn ~ 8%, where isolated Mn^+ spins com- 
prise ~ 2.9% of all cation sites. In Fig. 4b we investigate 
whether it is then possible - with the same total num- 
ber of Mn^^ spins - to increase the number of isolated 
spins by redistributing the Mn^+ in digital planes (solid 
dots) . Bulk can be thought of as 2D planes of spins with 
x\f.^=8%, spaced every monolayer. Next, we consider 
2D planes with twice the density {xj^^ =16%), spaced 
every other monolayer, which results in a paramagnetic 
enhancement of over a third, as shown. However, spacing 
planes with a;|P^=24% every third monolayer results in 
fewer free Mn^+ spins per unit volume than in the case of 
bulk. Additional divisions continue to reduce the para- 
magnetic response. So, only by spacing magnetic planes 
every other monolayer is it possible to increase the den- 
sity of free Mn^+ beyond 3D spin distributions. However, 
any intermixing during growth couples the 2D planes and 
dramatically reduces the paramagnetic enhancement, as 
shown by the open dots for the case of full segregation. Of 
course, clever schemes for control of the spin distribution 
within the 2D plane could certainly result in reduced spin 
clustering, such as MBE growth in the (120) direction, 
where neighboring-Cjation sites in the (120) plane are not 
nearest neighbors.E3 Thusfar, however, such efforts have 
been hampered by the inevitable inter-monolayer mixing 
of atoms during growth, leading to spin clusters. 

We apply the model to measurements of superlat- 
tices and quantum wells containing "digital" planes of 
DMS. Fig. 5a shows the measured splitting between ex- 
citon spin states in two superlattices with nominally sin- 
gle monolayers of Zn.75Mn.25Se and Zn.50Mn.50Se (sep- 
arated by 4ml of ZnSe). The dotted lines are Bril- 
louin fits to the low- field magnetization {H < 8T). In- 
creased spin clustering in the Zn.50Mn.50Se monolay- 
ers is evident in the smaller paramagnetic saturation, 
and more linear high-field susceptibility. With perfect 
2D planes, however, it is impossible to account for the 
15% larger paramagnetic saturation from the superlat- 
tice with Zn.75Mn.25Se planes. However, assuming expo- 
nential, segregated Mn^+ profiles oc e^"/'^ for each of the 
Zni_a;Mn3;Se planes (reasonable for the low growth tem- 
perature of 300 C), the relative low field saturations can 
be reproduced with a decay length A= 1.1 5ml, implying 
partial segregation during growth. As a final study (Fig. 
5b), we attempt to account for the size of the magnetiza- 
tion steps observed in PL from the quantum well contain- 
ing twelve l/4ml planes of MnSe. Magnetization steps 
arise from the partial unlocking of antiferromagnetically- 
bound Mn-Mn pairs, resulting in a step height propor- 
tional to the number of pairs. The observed magne- 
tization steps are never more than 5% of the low-field 
'saturation magnetization' Mgat, a ratio which is much 
smaller than predicted by any conceivable distribution 
profile of the Mn^+ within the quantum well. The ex- 
pected step height for 3D, 2D, and segregated 2D spin 
distributions are shown for comparison. This puzzling 




Field (Teda) Field (Tesla) 

FIG. 5. a) Measured spin splitting at 1.6K for the two su- 
perlattices shown, from reflectivity data. Dotted lines are 
Brillouin function fits to the low-field data, b) The measured 
magnetization steps (via PL) in the quantum well with 1/4 
ml magnetic planes. The steps are smaller than any Mn^"*" 
distribution would predict. 

anomaly is seen in all 'digital' samples, and even quan- 
tum wells containing bulk {x^}^,^ =8%) DMS show a sim- 
ilar deficit. We postulate this effect is due to the nature 
of the PL measurement itself, which is not a direct mea- 
sure of magnetization, but is rather only proportional to 
the magnetization through the Jsp-d exchange interac- 
tion and the Mn^+-exciton wavefunction overlap. It is 
anticipated that true magnetization studies will reveal 
the correct magnitude of the magnetization step. 

In summary, we have presented a method for calculat- 
ing the exact number of spin singles, pairs, triples, and 
higher order clusters for an arbitrary magnetic concen- 
tration profile in the (100) growth direction, to model the 
magnetic properties of real, quasi-2D spin distributions 
in DMS heterostructures. Calculation of the magnetiza- 
tion for diffusion and segregation profiles reveals nearly 
identical values, so that fitting an intermixing length is 
likely possible only when the form (exponential, gaus- 
sian, etc) of the quasi-2D profile is assumed a priori^ 
as was demonstrated for the case of ZnMnSe:ZnSe su- 
perlattices. The model also predicts a larger paramag- 
netism compared with bulk spin distributions only if digi- 
tal planes are spaced every other monolayer, although the 
effects of intermixing will greatly reduce any enhance- 
ment. Lastly, the discrepancy between the magnitude 
of observed and predicted magnetization steps remains 
outstanding. The methods outlined in this paper will 
be of use in modeling future epitaxially-grown DMS het- 
erostructures, where spin distributions can be engineered 
with nearly monolayer precision. 
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